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Abstract 

. A study of the relationships between large-scale cloud fields and large-scale circulation patterns is 
presented. The basic tool is a multi-level numerical model comprising conservation equations for tem- 
perature, water vapor and cloud water and appropriate parameterizations for evaporation, condensation, 
precipitation and radiative feedbacks. Incorporating an equation for cloud water in a large-scale model is 
somewhat novel and allows the formation and advection of clouds to be treated explicitly. The model 
is run on a two-dimensional, vertical-horizontal grid with constant winds. It is shown that cloud cover 
increases with decreased eddy vertical velocity, decreased horizontal advection, decreased atmospheric 
temperature, increased surface temperature, and decreased precipitation efficiency. The cloud field is 
found to be well correlated with the relative humidity field except at the highest levels. When radiative 
feedbacks are incorporated and the temperature increased by increasing C0 2 content, cloud amounts 
decrease at upper-levels or equivalently cloud top height falls. This reduces the temperature response, 
especially at upper levels, compared with an experiment in which cloud cover is fixed. 
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The compensating features of the cloud greenhouse and the cloud albedo effects, allied to possible 
changes in cloud height, distribution, and characteristics complicate our understanding of the response 
of the climate to external stimuli, such ns increased CO 2 levels, The problem Is twofold. First it is 
necessary to understand the radiative effects of cloud cover, This has been discussed by Ohring and 
Clapp (1980), Cess and Rnmanathan (1978), Schneider (1972) and many others, The second, and less 
well understood, problem lies in understanding tho causes of cloud variability, or alternatively predict- 
ing cloud response in a given situation. Predicting cloud response docs not simply involve predicting 
total cloud cover. For example, Schneider has demonstrated the importance of the cloud top height. If 
cloud top height increases then the effective emitting temperature decreases and the greenhouse effect 
is enhanced. On the other hand, if total cloud area is simply increased then it is likely that the albedo 
effect will outweigh the greenhouse effect (e.g,, Hartmann and Short, 1980), However, Cess (1976) 
and Cess et al (1982) argue that the competing cloud effects may compensate each other very closely at 
all latitudes and seasons. Presumably, this docs not mean that effects at individual longitudes will also 
be closely compensated (Shukla and Sud, 1981). 

Much has been written about the physical characteristics affecting cloud formation and associated 

precipitation (for a comprehensive account sec Ludlum, 1980). The task facing climatologists is not 

only to understand where and why clouds of various shapes, sizes and characteristics form but to 

parameterize this in terms of the large scale. Although simple parameterizations of cloud cover, such 

as setting it a function of relative humidity (first suggested by Smagorin'sky, 1960; see also Slingo, 

1980) or allowing cloud where condensation occurs, appear at first glance to give climatologically 

acceptable results in GCM’s, high clouds are poorly treated and hence even qualitative uncertainty in 

the radiative feedbacks must arise. Furthermore, it seems that the parameterizations must be tuned to 

tho particular GCM used, clearly an unsatisfactory procedure. 

♦ 

At the risk of oversimplifying cloud physics, let us assume that clouds occur primarily where par- 
cels of air become supersaturated with respect to a plane surface. If the rate at which supersaturation 
occurs is sufficiently large to produce sufficiently large amounts of cloud water, then coalescence and 
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accretion of various cloud drops produce still larger drops which eventually fall to the surface. Tran- 
sport processes will extend clouds into regions (both space and time) where cloudwater evaporates. 
While the above notions are partially understood in specific cases, the associated large-scale processes 
that will produce regions of supersaturation and subsaturation are numerous and not well understood. 

A number of workers, including Manabe and Wetherald(1975), Roads(1978a), Schneider et o/(1978), 
and Wetherald and Manabe (1980) have shown that middle tropospheric clouds tend to decrease with 
increasing temperature even though the intensity of the hydrologic cycle (precipitation and evaporation) 
increased. The explanation is that with a more intense hydrologic cycle, the vertical velocity variance 

- tends to increase, resulting in more precipitation in regions of ascent and more drying in regions of des- 

» * ‘ 

cent. In regions of descent where cloud cover was already present or under conditions where horizontal 
advection of dry air is important, the drying effects have more effect on the cloud cover than the 
moistening- effects and hence cloudiness decreases with increasing temperature. 

There are certainly other factors influencing the cloud fields. Roads (1978b), for example, has 
shown that under certain conditions zonal mean cloud cover can be decreased by increasing the mean 
downward zonal velocity, by increasing the diabatic heating, by increasing the static stability and by 
decreasing the scaled evaporation, q y . Scaled evaporation is the ratio of the surface saturated humidity, 
ft', to an atmosphere saturated humidity, q' a . If temperature alone increases, but lapse rate remains the 
same, the scale evaporation decreases, since B/dT(q y ) <0 for T s > T a . 

In a GCM experiment, Wetherald and Manabe (1980) found that low level fractional cloud cover- 
age actually increased with increasing temperature in high latitudes. (Similar results were also found by 
Potter et al (1981)), This occurred over a region in which surface albedo decreased rapidly, due to the 
poleward retreat of the surface ice sheet, causing surface temperatures to rise noticeably - a rise of 19K 
was found at latitude 80° in Wetherald and Manabe’s model due to a 6% rise is solar constant. Atmos- 
pheric temperatures did not rise as rapidly, because of the moderating influences of cross-latitude heat 
exchange'. Vallis (1982), Roads and Vallis (1982) explained this using relatively simple zonally aver- 
aged models. As the ice sheet retreated, the lower atmosphere (in high latitudes) was destabilized and 
an increase in scaled surface evaporation ensued which caused cloud cover to increase. In lower lati- 
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tudcs however, where surface albedo stayed fixed, Roads and Vallis found that the lapse rate was stabil- 
ized and sealed evaporation fell. Acting alone this would cause cloud cover to fall; however this was 
compensated for by an increase in the net radiative cooling required to balance the increased precipita- 
tion. 

Several mechanisms have thus been identified that affect the large-scale cloud and relative humi- 
dity fields: (1) eddy vertical velocity, (2) horizontal advcction, (3) low-level temperature inversion, (4) 
radiative effects. How important these mechanisms are for models with more verisimilitude to nature is 
still unknown. For example, one of the chief simplifications (excepting the GCM work) has been the 
restriction of clouds to one-level, so prohibiting understanding of the vertical structure of the cloud and 
relative humidity fields. Since the vertical extent of cloud fields is crucial in determining the relative 
importance of the solar and infrared radiation components of the radiation field this simplification is 
ultimately unacceptable. Nor docs any wholly satisfactory scheme currently exist, relating cloud cover 
and cloud optical properties to the dynamical fields (c.g., relative humidity fields) predicted by a GCM. 

The study to be described in this paper addresses the problem of the vertical cloud structure using 
a two-dimensional model with variations in the zonal and vertical directions, The problem of relating 
clouds to the dynamical and thermodynamical fields is addressed by including an explicit equation for 
cloud water. The model is written in such a manner that the condensation and radiation algorithms 
could easily bo added to a more comprehensive general circulation model. Some of the suggested 
hypotheses regarding cloud variability could then be tested further. 

2. Model 

The model comprises conservation equations for temperature, water vapor, cloud water and 
paramctcrizntions for the appropriate sources and sinks. The continuous form of the equations may be 
written; 

Thermodynamic 

It + v -° + " «^ F » + 1 -kc«y~* + xv'o 


a) 


■ - 6 - 


Watcr vapor 

^2. + Vv<7 + ~Q^ (,> q ** 8~^Fq + C q + K'7 2 q (2) 

Goud water 

4? + V*vc + -l-toc - ^J-F c + C c + KV 2 c . 0 } 

dt ~ dp dp 

Precipitation (steady state form) 



where 

0 «= j potential temperature = Tp~ K 

i 
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v *=» horizontal velocities 
<y «= pressure velocity 
g *= gravity 

C p *= heat capacity at constant pressure 
L = latent heat of condensation and sublimation 
q *= water vapor mixing ratio 
c «= cloud water mixing ratio 
t *= time 

V* — horizontal divergence operator 
p *= pressure 

F g = thermodynamic eddy + radiative fluxes 
F q = water vapor eddy fluxes 
F c *= cloud water eddy fluxes 

r * 

F/ «= precipitation fluxes 
C *= conversion terms for water substance 
K *= horizontal diffusion coefficient = 10 4 /n 2 s“ l 
R «= perfect gas constant 


k - RfC p 
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2.1. Numerics 

The continuous equations are finite differenced in the horizontal and vertical. The vertical grid 
consists of points (9) equi-distant in pressure between 1000 and 100 mb. Fictitious boundary points 
below 1000 and above 100 mb, are used to define values (via boundary conditions) at 1000 and 100 
mb. These boundary conditions at 100 mb are 


or 

Op 

0p 


0 

0 


0c 

Op 


«o 


0 ) 0 

and at 1000 mb 


i 


k^-k s (e s -e) . 

- pk s (.q B -q) 

k~- **~k s c 
Op 

6)“0 

where /3, k and k s are empirical coefficients and will be explained later. 


Sjxtecn equally spaced gridpoints are used to represent variations in the zonal direction, and 
periodic boundary conditions are imposed. In the north-south direction the model currently assumes 
no variation. 


The equations are stepped forward in time using a leapfrog and periodic restart scheme for the 
dynamical terms and a forward step for the parameterizations. It was found necessary to filter the tem- 
perature and moisture in the horizontal to remove 2 gridpoint waves. The condensation routines are 
then called in orJer to maintain convective stability and remove any supersaturation. 

The velocity fields are preset and remain constant throughout the integration. In addition to a con- 
stant zonal wind of constant shear, a divergent wind is also included such that with the pressure velo- 
city the mass continuity equation is satisfied That is 


* 8 


Vv-- 


dco 

dp 
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atjd w-0 at the top and bottom is imposed. The imposition of two-dimensionalily and fixed velocity 
fields are the primary simplifications of this model. The paper thus becomes a study of the relationship 
of cloud structure to a given dynamical field. Complete verisimilitude is sacrificed in the interests of 
isolating perticular mechanisms. Clearly the dynamical effects and feedbacks of cloud systems on the 
large-scale must be neglected (e.g., Moncrieff and Green, 1972; Stevens and Lindzen, 1978). 

2.2. Vertical eddy fluxes and diabatic heating 

The surface fluxes of heat, moisture and cloud-water are described by simple bulk aerodynamic 
laws . 

! 

F e - k s (0 g —0 ) 

F q - 13k, ( qg-q ) 

I'c » k, ( — c) 

where /3, which h set to unity here, is related to soil moisture. Further up in the atmosphere the 
diffusive fluxes are • , 

i 

F 0 -k^- 

dp 


F e - k 


A 

dp 

9c 

dp 


The boundary condition used at the surface is that these two fluxes are equal. 

The diabatic heating due to radiation is initially parameterized by a Newtonian cooling law 

g^-Fg(rad) - (0 e -0)/lO 6 s 
OP 

More elaborate radiative parameterizations (i.e., Lads and Hansen (1976) for solar and RodgerS (1967) 
for infrared) will be incorporated in Section 4. 


2.3. Water conversions and convection 

The cloud water conversion and precipitation parameterizations arc essentially based on 
simplifications of physically based schemes. A comprehensive treatment of the more complex theory is 
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given by Kessler (1969). Here we assume that elaborate nvtcrophysical and eddy transport schemes arc 
unnecessary and that the transports of water vapor and temperature are accomplished by large-scale 
advection or can be adequately parameterized by a convective adjustment scheme. 

The conversion terms, C, are complex nonlinear terms that describe the conversion between 
water vapor and liquid water or ice and the conversions between different number densities and size 
distributions, The conversion in this model is simplified by the perhaps rather arbitrary distinction 
between watcr/ice droplets of very small diameter and negligible fall velocity (cloud water) and 
watcr/ice droplets of large diameter and a substantial fall velocity (precipitation). 

2.3.1. Conversion of cloud water to precipitation 

The following processes are thought to be the dominant contributors to the conversion of cloud 
water to precipitation. Initially, small cloud droplets self-collect to form larger cloud droplets and small 
precipitation particles in a process referred to as auto-conversion by Kessler. The process seems to 
have a cubic dependence upon cloud water. The so-called Bergeron process also contributes to the for- 
mation of small precipitation particles by augmenting ice particles at the expense of liquid water in addi- 
tion to sublimating water vapor. Presumably this process decreases the auto-conversion time scales. 
Once small precipitation particles are formed, they fall through the cloud and further collect cloud dro- 
plets. This is referred to as accretion and is dependent upon the amount of cloud water and precipitable 
water as well as their respective number densities and size distributions. In the limit of large accretion, 
cloud water is probably converted to precipitation as the 4.5 power of cloud water. 

It would thus seem that precipitation is converted to cloud water in a highly nonlinear fashion 
(see also Sundqvist, 1981). This insures that clouds have a fairly uniform mixing ratio (compared to, 
say, the mixing ratio for water vapor) that precipitation falls from only dense clouds, and that drastic 
changes in the circulation parameters are only likely to give small changes in the cloud water. For this 
study we feel it is adequate to use a simple conversion rate, namely 

C ct “ ac 3 (5) 

where a “ 10 4 s _1 . If the precipitation rate is typically 1.9.vlO“ 5 A'g/;r 2 s _I then the average cloud water 

mixing ratio is 5.7xl0~ 5 kg/kg. If clouds exist in only 50?/o of the atmosphere then the average cloud 
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water content in the clouds is 7.2xlO -5 and the average cloud water content (cloudy ■+■ noncloudy 
regions) is 3.6xl0“ s . These estimates of cloud water content are all within Sasamori’s (1975) esti- 
mates. The multiplying factor and exponent will be tested later to determine the sensitivity of the 
model to the cloud physics parameterization. As will hs shown, the parameterization is crucial for 
determining the characteristics of the cloudwater but negligible for determining the cloud fractional 
cover and large-scale vapor fields. This occurs chiefly at the lower levels because the amount of cloud- 
water in the atmosphere is negligible compared to the amount of water vapor. 

Because of the highly nonlinear rate at which cloud water is converted to precipitation and the 
desire to maintain a fairly large time-step sufficient only for large-scale dynamical stability, it is assumed 
that the process can be described as 1 

f+A/ 

C c j ■=» — J ac l dt 

• t 

or 

Cci “ ~r~ (1 — j- — r~ — ^r; t) (6) 

c! At (1+2A tac 2 )* . • 

2.3.2. Conversion of precipitation and cloud water to water vapor 

Once precipitation is formed it falls through the cloud and into the, free atmosphere. If the pre- 
cipitation is falling in a region of unsaturated air then the precipitation is assumed to evaporate com- 
pletely or to the point where the air is just saturated. That is, the following set of equations are solved 


C P ST+ Lhq -0 

q -f Sq «=» q ( T,p ) + 

Ap 

where 8 denotes the differential changes and an 


denotes saturation values (saturation is over water for 


temperatures above freezing and is over ice for temperatures below freezing). The conversion of pre- 

« 

cipitation to water vapor is 


- , (< 7 ‘-<?) 

-k r s • 

A ' (i+ #fr> 
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If this results in — (/} 4* 8F/) < 0 then the precipitation is completely evaporated without saturating the 
atmosphere and the conversion is 


t Fl ■ 

Cloud water in unsaturated conditions is treated in a similar manner to give C C q» 


2.3.3. Conversion of water vapor to cloud water 

The conversion of water vapor to cloud water requires simply that supersaturation is reduced to a 
saturation value instantaneously; hence for supersaturation conditions the following set of equations is 

i 

. solved 

CpST 4* 0 

q + l4«*q m (T,p) + &?8T 

or 

c k=al 

v C p BT ... ( 

2.3.4. Convective condensation 

Ostensibly the conversion terms have now been parameterized for the model, but one process 
remains that will alter the condensation process, namely moist convection. As stated by Manabe et a! 
(1965) and Sarachik (1980) our knowledge of this process is insufficient to warrant a too detailed 
parameterization of moist convection. Further, Lindzen et al (1982) found that sensitivity experiments 
with a moist-adiabatic adjustment yielded similar results to experiments with an alternative cumulus 
cloud model. Thus, we employ a moist convection adjustment scheme, which operates when the verti- 
cal gradient of moist static energy is negative. In the adjustment, it is assumed that turbulent processes 
occur which equilibrate the moist static energy between the levels, maintain saturated conditions, con- 
serve the total energy, and equilibrate the cloud water. Momentum could also be equilibrated in a simi- 
lar fashion. The equations for two contiguous levels (1 denotes the upper level and 2 the lower level) 


are: 
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c p bt 2 + q,8 r, + L5(?2 + - o 

A/r + -%(pr + - •“} + a ( 8(72 ~ 5<7t) *" o 

2 pf pf 

$2 + 8tf2 ” 92 + "ay 5 ^2 

<7l + 8(71 ** <7i + *^y , 87’i 
8cj + 8C2 “ + 8^) 

d + 8cj “ c 2 + 8c 2 

where the difference in moist static energy between two levels has been defined as 
A/r - -^(/>f + Pi) (02 - 0i) + Uq 2 - <7i) 

The convective process occurs whenever Ah > 0 and (ci 4- 5c\ + c 2 + 8C2) ^ 0, The process in the 
model begins by testing and solving contiguous levels from the bottom of the atmosphere and moving 
upward. Although more than two levels could be solved simultaneously, idealized tests showed that 
iteration of the above two-level process produced the same results as multi-level systems; hence, for 
simplicity, only the two-level system is used. Also, since the time stepping was sufficiently small 
(approximately 1 hour), iterations at each time step produced the same results as no iterations. How- 
ever, if the condensation routines are called infrequently, it becomes necessary to iterate. For rapidly 
changing dynamical conditions it may also be necessary to iterate. 

To summarize the above processes, then, the model initially tests for moist static instability. If 
instability exists the temperature and moisture content of two levels are adjusted to equalize the moist 
static energy and convert excess water vapor to cloud water. The convection is assumed to equalize 
cloud water mixing ratio in the vertical. This is repeated through the depth of the atmosphere. Cloud 
water is then advected, and converted to precipitation according as (6). The precipitation then falls, 
either to reach the surface or to be evaporated in non-saturated layers. If ,10 convection occurs, we 
may still have local cloud water vapor to cloud water conversion, and the precipitation process follows 
in precisely the same fashion. 

Having calculated the new stable state, it becomes necessary to define the condensation and the 
accompanying turbulent transports. These are given by 
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dF n sr, 

HF a/ 

^fl/r2 9^2 La 

*17 " IT ~ -c, c * 

S/fyi 0<7s 

^1T" a/ 

„ $3l + c 
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where 


A (8<7rH5tf 2 ) 

. * " ~ Ar 

Q 

Dry convection also occurs infrequently whenever As *■» -y'(p{ : +pf) (0 2 -0i) > 0. As in the 

moist convective process, q and c are assumed to be well mixed while conserving the dry static energy. 

Havn>§ found the various functions, ‘the equations (1, 2, 3) for the time rates of changes due to 
th-4 conversion and convective terms and the precipitation equation (4) become 

If- - sfjf„ + - Q,-C„))r* 

If “ *17 ~ & + + 

0c JFc i A A _ A 

dt™ 8 dp + Cflc Ccq Cc/ 

-/}«/ ( C c ,-C lq )(ip/g 


3. Experiments (Newtonian cooling) 


The experiments with a simple Newtonian cooling scheme are designed to illustrate the effects of 
changing the horizontal advection, vertical velocity and surface temperature, and to study the effects of 
incorporating an explicit cloud water equation. The velocity holds consist of a mean zonal wind, plus a 
zonaMy varying vertical velocity wave. The velocity fields are not meant specifically to reproduce a par- 
ticular dynamical situation, but to aliosv. the effects of horizontal and vertical motion to be studied in a 
simple, framework. 


la* 


Figure 1 shows the vertical profiles of the various constants for the control run. The vertical 
velocity is a maximum (1,6x10“*^ m"' s~ 2) at 550 mb and decreases to zero at 1000 and 100 mb. 

The vertical mixing parameter k decreases from 12.6xl0“ l at 1000 mb to 16.5xlO -3 (kg 2 m ~ 3 s~ 3 ) at 
100 mb. The zonally averaged wind increases from 10 at 1000 mb to 28 (ms~ x ) at 100 mb. Since the 
divergence is assumed to be zero, the zonal wind also includes a component that varies in x. This Com- 
ponent is zero at 550 mb and increases linearly to a wave with an amplitude of 4.46ww“ l at 1000 and 
100 mb. If the pressure velocity varies as a cosine wave, then u varies as a sine wave with positive 
amplitude below 550 mb and negative amplitude above. 

Because of the constant value of the winds, steady state solutions for the equations are obtained. 

* 

We will describe in detail the results for the control experiment and then results from other experi- 
ments motivated in part by the simple one-level cloud system of Roads (1978b). 

Figure 2 shows a steady-state solution as a function of longitude for the control parameters. Here 
the periodic domain extends from west to east with a distance of 5xl0 3 km and grid points every 312.5 
km. Essentially wave #6 phenomenon is studied if these scales are appropriate for midlatitude circula- 
tion systems. The vertical velocity wave is shown in the upper part of the figure to have maximum 
descending motion on the edges of the domain and maximum ascending motion in the center of the 
domain. The precipitation at the surface is correlated with the vertical motion field although precipita- 
tion is certainly narrower, occuring significantly over less than 1/3' of the domain. Note also that the 
precipitation occurs almost strictly on the leeward side of the vertical velocity wave. The characteristic 
features suggest that for the chosen parameters, horizontal advection of moist air is extremely impor- 
tant for the precipitation field. 

Various measures of the total cloud cover are also shown in Figure 2. Of course, the most natural 
criterion to use is simply the presence of cloud water, although clearly this cannot be used in GCMs 
which do not predict cloud water. The narrowest measure is the area of positive condensation which 
also corresponds with the area of precipitation at the surface. If the presence of cloud water or a relative 
humidity criterion (H=T) is used then the cloud boundary extends a further distance on the windward 
side. Smaller relative humidity criterions naturatly give larger coverages still. 
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The two-dimensional structure of the cloud, relative humidity and temperature fields are shown in 
Figure 3a,t>,c. First note that the cloud has an eastward tilt with altitude indicating that dry air advec- 
tion occurs on the windward side and cloud water advectlon occurs on the leeward side. The relative 
humidity field gives a good indication of this structure and also shows that a minimum in relative humi- 
dity exists below the tropopause. The temperature field is that which would be expected for a forced 
circulation with the strongest positive deviations eastward of regions of descent and the strongest nega- 
tive deviations in regions eastward of the ascending regions. The deviations are substantial with a max- 
imum eastward temperature difference of about 15°C indicating differences in heating of 1.5^10“ 5 ^fs -1 . 
Of interest also is that the deviations are negative in the region of the cloud top indicating that the 
parameterized heating rate is positive in this region. 

. Various experiments were run to see the contributions and sensitivity of various mechanisms. 
Figure 4 shows the response when the mean zonal wind is removed. In this case tb** cloud cover 
decreases at the lowest level and substantially increases at the upper level. The relative humidity 
becomes lens in the descending regions and the atmosphere becomes relatively warmer in the descend- 
ing regions. Not much temperature change occurs in regions of ascent except at the top of the cloud 
where strong negative temperature deviations occur indicating strong radiative heating. Due to this 
drastic change in circulation pattern, the change in the zonally averaged relative humidity docs not 
correspond to the change in cloud coverage. For example, in the middle levels fractional cloud cover 
increased while the zonally averaged relative humidity decreased. 

Another experiment in which the vertical velocity amplitude was decreased an order of magnitude 
yielded cloud coverage at all longitudes and levels except for the uppermost level and part of the next 
uppermost level. Although the decrease in vertical velocity yielded eloud eoverage and evaporation 
everywhere, the precipitation is much less for this case as may be seen in Figure 5. Note also that the 
precipitation is strongest. in the center of the rising motion and much weaker in regions of descent, 
decreasing to zero in regions of strong descent. Thus the cloud and precipitation characteristics are 
strongly dependent upon the ability of the atmosphere to dry itself out in addition to being able to pro- 


duce clouds. 
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Another mechanism that is important for latitudinal differences in cloudiness is the mean zonal 
vertical velocity. Although its inclusion is inconsistent because meridional divergence is not allowed, a 


term of the form was added to the flux V-v/t + -—a)/!, where A is temperature or mixing ratio. 

dp ~ dp 

5 had the same vertical profile as the eddy vertical velocity term (see Fig. 1) but a much smaller ampli- 
tude maximum of 2.5xl0“ 2 kgm~'s~*. Because the continuity equation is no longer satisfied, precipita- 
tion was much less than evaporation. However, a steady state was reached and is shown in Fig. 6. 

Note that a low-level cloud of small areal coverage is formed in a very dry and hot atmosphere. 


When the surface temperature is varied, large changes in the cloud structure result. Fig. 7a, b 
show the response in the cloud and relative humidity fields for an increase of surface temperature of 
10K and Fig. 7c, d show the response in the cloud and relative humidity fields when the surface tem- 
perature is decreased by 10K. Note that clouds strongly increase for a more unstable boundary layer 
almost everywhere. Also of interest was that the clouds were almost completely formed by moist con- 
vection for the unstable boundary case and completely by local supersaturation condensation for the 
stable boundary case. As discussed by Roads and Vallis (1982) these cases have applications to models 
in which the ice line is allowed to change, because of large changes in low-level static stability and 
cloudiness in regions of the ice line. 


It was pointed out above that a major difference between this model and others is the explicit 
inclusion of cloud water. If all cloud water is precipitated out as soon as it forms, then essentially the 
model reduces to the conventional types. A comparison of cases with and without cloud water and its 
concomitant advection is useful. Thus three of the above described experiments were rerun with 
immediate precipitation of cloudwater. For these cases it is shown in Figure 8 that the relative humi- 
dity fields are more or less the same with or without cloud water except at the top of the atmosphere 
where cirrus cloudiness is dependent upon cloud water being present. Because of the large influences 
that cirrus cloudiness has on the radiation field it would seem to be better to retain the cloudwater 
parameterizations. However, good estimates of cloud cover in the lower regions are obtained by equat- 
ing cloud cover to values of relative humdity fields near unity. At the other extreme Figure 9 shows 
the response when C c i « i0 2 c 3 (as opposed to lO^c 3 see eq. 5). The cloud wafer content increases by a 
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factor of 5 with the largest areal change occurring at the upper levels where cloud water can now extend 
over much more extensive regions because of the importance of advection. The change in the relative 
humidity and temperature field are negligible however at the lower levels. Thus, from the two previous 
experiments, we conclude that cloud characteristics and variations are determined by the circulation 
parameters that result in saturation and that modifications to these occur chiefly at upper levels where 
cirrus clouds detrain from the main cloud. The precise form of the cloud microphysical constants 
appears unimportant. 

The change in the Newtonian heating profile produced some interesting effects. Figure 10 shows 
the response when the equilibrium temperature was either decreased by 10° everywhere, including the 
surface, or increased by 10°, including the surface, everywhere. Note that the cloud fields remain 
about the same in the lower to upper troposphere but that the run with the increased temperature 
resulted in decreased upper-level cloudiness. 

Somewhat different responses occur when the cloud water is totally removed at each time step, 

Fig. 10a and 10b, because then the amount of upper level cloudiness (//> 1) is actually decreased at 
the upper levels (except for the topmost level) with decreasing temperature. So if cloudwater is 
included, the upper levels tend to detrain cloudwater over a longer distance. Presumably this occurs 
because at the higher temperatures the relative humidity deficit is sufficient to evaporate the cloud 
water but at lower temperatures the cloudwater (which is maintained at a fairly uniform level by the 
nonlinear precipitation processes) is much larger than the relative humidity deficit and can travel over 
much larger distances before finally being evaporated in the driest regions. Thus, if cloudwater is 
retained in the model, upper level cloudiness will increase with decreasing temperature whereas if cloud 
water is not retained, upper level cloudiness will increase with increasing temperature. 

4. Cloud Radiation Feedbacks 

Schemes for calculating the infrared and solar radiation fluxes were incorporated in order to allow 
a more self-consistent model. The schemes are fully interactive in that the water vapor and cloud water 
content predicted by the model are used in the radiative calculations. 


-is- 


4.1. Infrared radiation 
* 

The infrared scheme is based on the emissivity approximation (see Rodgers, 1967), in which the 
integration over all spectral wavelengths is incorporated into simple transmission formulas. In terms of 
transmissivity the equations are 

L 

F\ (p) “ J J9(pVr (/>,//) 4- B s T(p s ,p) 

r(p t p) 
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r(p.O) 

B(p')c1t(p,p') 

1 

where 


F\(p) 

== the upward flux 

Flip) 

*= the downward flux 

r (/>,/>') 

= transmissivity between 

B(p ) 

*= crT*(p ) 

Tip) 

— absolute temperature 

C r *=* 

Stefan-Boltzmann constant 

B g 

- <r Tg 


T g = surface temperature 

The transmissivity for water vapor is taken from Rodgers. Carbon dioxide is accounted for by assum- 
ing that the total transmissivity is the product of the transmissivity of H 2 0 and C0 2 where 

r(CC?2) “ .926 —.02 /«(Am) 

Am - 120^- 
Ps 

where p s is the surface pressure and A p 7 is the pressure squared differential (using a value of 120 
implies a C0 2 content of approximately 300 ppm). Thus Am is a pressure weighted path length. In a 
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similar manner, the transmissivity in the presence of cloud is the product of cloud, Hi 0 and CO% 
transmissivities. Currently the cloud transmissivity is taken to be 

r «" I u/a\u<a 

t *" 0; 

where 

dti ** 

8 

a *• 3.vl0*’ 2 A v e nr 7 

Tiro dependence of r on u produces nearly black clouds everywhere except at the cirrus levels. For 
. example, at the uppermost level the typical c is on the order of LvtO ~ s kg/kg and the typical r is on the 
order 2/3, t 

4.2, Solar radiation. 

t 

The solar radiation scheme uses the multiple scattering method described by Lacis and Hansen 
(1974) and Somerville ct al (1974). Simplified two stream approximations and absorption coefficient 
approximated by a sum of C> gray absorbing intervals arc used in partially cloudy skies. In an entirely 
clear sky and above the highest clouds the absorption is calculated using an integrated absorption 
coefficient. In this study single layer optical depths are assumed to be given by 

t *■ e*”^~ 10 2 
8 

The typical optical depth in the lower layers is then around 8 for a pressure interval of 100 mb and a 
cloudwatcr mixing ratio of 8\iO“ s kg/ kg. In the upper-most layers the optical depth is about 1 for a 
pressure interval of 100 mb and a cloudwatcr mixing ratio of LvlQ~ 5 kg/kg> 

As a typical cosine zenith angle we use /a ** .225 which is appropriate for the annually and diur- 
nnlly averaged zenith angle at 45° latitude. The solar constant is set at 1372 IKw 2 . 

4.3. Surface temperature 

» 

The surface temperature equation can be written in the form 



4 . 4 . Experiments j 

I 

t 

The equilibrium profiles of cloud water, relative humidity and temperature are given in Figure 11. j 

i 

The mqjor differences between the previous cases is that the lower boundary tends to be more unstable 

4 

* 4 ? 

which produces cloud at the lower levels everywhere the motion is upward. The upper level cloud field 

is similar however. • * ; 

! 

, i 

^ 0 , ' | 

The radiation cooling rates, — , are shown next in Figure 12. Immediately noticeable is that in ! 

i 

the atmosphere infrared radiation is dominant especially near the cloud boundaries. Above the lower 

i 

level clouds, cooling rates of up to 157day can occur whereas below the upper level clouds heating t 

i 

rates up to 6°/day can occur. In contrast is the surface flux balance shown in Figure 13. Here it is the j 

| 

solar flux that varies most strongly from east to west with substantially less radiation below the clouds. j 

This solar radiation variation is balanced by an increase in downward radiation, a decrease in upward j 

i 

radiation and decreased sensible and latent heat fluxes to the atmosphere. 

Figure 14 shows the response when the cloud radiation feedbacks are modified. Here the model 
was only integrated for 30 days and although the solution had not yet reached a steady state the solu- 
tion was in a statistically stationary state and was converging slowly toward equilibrium with a slowly 
increasing exponential damping rate of >30 days. Note that in the absence of cloud radiation feed- 


r ? 1 




1 


ft. 
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backs the temperature increases almost everywhere except at the upper levels. In the presence of 
infrared radiation feedbacks alone the atmosphere is still warmer, demonstrating the cloud greenhouse 
effect. Solar radiation, however, dominates the total cloud radiation feedback and this model atmo- 
sphere tends to be colder because of clouds. These results are similar to results found with the GLAS 
GCM (Herman et al, 1980). 

Finally, we show a calculation for the response of the model atmosphere to a quadrupling of CO 
This amourd of CO2 increase is well within current projections (Bolin et al, 1981) which show that by 
tne end of the ?,3rd century the atmosphere could contain 4 to 8 times present amounts. Two expsri- 

» » 

ments were run to test the effect of the clouds on the solution. In the fixed cloud, FC, experiment, the 
clouds’which interacted with the radiation field were the clouds found from the control run. In the 
/ - ' variable cloud experiment the clouds were allowed to change and interact with the radiation field. The 

time variation of the two runs for an arbitrary surface temperature is shown in Figure 15. Note that 
the variable cloud experiment produces a shorter convergence time than the one with fixed clouds, 
indicating a negative feedback. 

* 1 

i The change in the cloud field is shown in Figure 1-6. Note that the upper level cloudiness tends to 

decrease, and cloud cover at the eastern and western side of the original system also tends to fall in 

I 1 

if spite of the general overall increase in cloud water in the center of the cloud. In fact, the one positive 

| * 

f region on the left edge of the cloud was substantial, increasing by 7 aT 0 -5 and strongly increasing the 

| • infrared cooling in this region. As shown in Figure 17, this change in cloudwater decreases the albedo 

of the atmosphere indicating that the quadrupling of C0 2 causes slightly more absorption of solar radia- 
tion for the variable cloud experiment than the fixed cloud experiment. However, because of the 
increased water vapor content the albedos are reduced in both cases. 

1 

Figure 18 shows the zonally averaged potential temperature response for a C0 2 quadrupling for 
the fixed and variable cloud experiments. In the variable cloud experiments the quadrupling produced 
a potential temperature change of about 2°K at the surface increasing to about 4K. at 50kPa and then 
decreasing in the upper troposphere with an overall increase in potential temperature. In the lower to 
mid troposphere the fixed cloud experiments showed similar albeit larger responses and in the upper 
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troposphere the responses are completely opposite with the fixed cloud experiments showing a very 
strong increase in potential temperature. 

In the fixed cloud experiments the vertical temperature structure tended to change according to 
changes in the moist adiabatic lapse rate. However, in the variable cloud experiments temperature 
responded similarly up to about 40kPa but then the temperature change decreased rapidly. Presumably, 
the decrease in upper level cloudiness has a strong influence on the results of the variable cloud experi- 
ments: reduction in upper level cloudiness allows radiation from lower in the atmosphere to more easily 
escape to space. 

i 

This experiment shows that assumptions about cloudiness can have a large influence upon the 
sensitivity of the upper troposphere, and to a smaller degree the lower tropospheric response. In this 
respect we are reminded of the experiments of Weare and Snell (1974) which also showed a reduced 
sensitivity to C0 2 because of- their cloud model assumptions. (See also Charlock, 1982; Paitridge, 

1980.) 

5. Conclusions 

This study has been concerned with the relationship of cloud fields with various features of 
atmospheric circulation patterns. In particular we have examined the effects of explicitly including 
cloud water and the effects on cloud cover of increasing temperature with increased atmospheric carbon 
dioxide. We have also considered the effects of variation? in horizontal and vertical advection. 

The effect of including cloud water is clearly seen only at the uppermost levels. At mid and lower 
levels, cloud cover is very well correlated with the relative humidity and condensation fields. Only at 
upper levels, where advection of cloud water over large distances can occur, does its explicit inclusion 
make a qualitative difference. Further, the cloud cover is not strongly dependent on the particular 
microphysical parameterizations used, although the cloud water content (and hence the albedo) of the 
atmosphere is. When the temperature is altered,* cloud cover generally responds in a similar fashion as 
relative humidity, except at upper levels. With cloudwater, the upper-level cloud cover tends to 
decrease with increased temperature because enhanced evaporation of cloud water prevents extensive 
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advcction, However, with no eloudwnter, relative humidity increased with increasing temperature at 
the upper levels. 

The incorporation of a fully self-consistent radiative model did not essentially alter these conclu- 
sions. The cloud field that developed was similar to the ones which used Newtonian cooling although it 
was slightly larger at the lower levels. It was shown that infrared radiation was most important for 
determining the atmospheric temperature structure whereas solar radiation was most important for 
determining surface variations, If the cloud radiation feedbacks were excluded then the model atmo- 
sphere warmed slightly albeit much less than if only the solar radiative feedbacks were excluded, This 
indicates that the infrared and solar cloud feedbacks were almost compensatory but that the solar radia- 
tion feedback was slightly stronger. 

On the other hand, in a sensitivity experiment in which the C0 2 was quadrupled the dominant 
feedback was the infrared component. Increasing C0 2 levels warmed the atmosphere and generally 
decreased upper level cloud cover, but increased cloud water content slightly at mid levels. Mid level 
areal cloud cover was little changed. The reduction in upper level cloud cover caused a reduction in the 
cloud-greenhouse effect (larger than any albedo effect) and reduced the total warming, compared with 
the case in which cloud cover was fixed. It would thus seem that the explicit incorporation of cloud 
water is necessary to predict correctly the response to C0 2 changes. In this model, clouds act to reduce 
the sensitivity of the model response to C0 2 perturbations. 

The relationship of the cloud cover to the dynamical fields are generally consistent with the previ- 
ous pictures obtained using one-level systems (e.g., Roads, 1978b). Thus, areal cloud cover is 
increased by decreasing eddy vertical velocity, increasing the mean upward vertical velocity, decreasing 
the zonal advcction. 

Finally we mention again some of the limitations of this model. The main simplifications lie in 
the assumption of two-dimensionality and in not allowing the large-scale cloud systems to change the 
winds in response to temperature changes. Thus the increased temperature - increased intensity of 
hydrology cycle - increased vertical velocity - reduced cloud cover feedback was not included in the 
CO 2 experiments. Further, the stratosphere was assumed to be an isothermal layer over an active 


troposphere. The relaxation of such assumptions within the framework of a simple model will be an 
important next step in understanding the cloud and climate response to changes in external forcing. 
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Vertical structure of various control parameters. A denotes a genera! com- 
ponent for either w, n-10, k y co and the multiplying factor is denoted on the 
figure. The temperature scale for T 6t 0 t is given at the top of the figure. All 

* 

values are in mks units. 

Zonal structure of various components as a function of I, the zonal gridpoint 
running from east to west, w denotes the eddy pressure velocity, P the precip- 
itation rate, and E the evaporation rate in kg m~ 7 s~ [ . Various measures of 
the total cloud cover; C, cloudwater > 0; C+, condensation > 0; 

//> 1; //>. 9 are also given. 

The two-dimensional structure of the cloudwater for the control case: Con- 
tours are every 0, 4, 8xl0 -5 kg/ kg. 

Same as 3a except for relative humidity. Contours cover the range 30% to 
90%. 

Same as 3a except for deviations of potential temperature field from the radia- 
tive equilibrium values. Contours are every 5 degrees and positive values indi- 
cate Newtonian cooling and negative values indicate Newtonian heating. 

Same as 3 except 77 = 0. 

Zona! structure of the precipitation and evaporation for the cases 77 = 0 and 
o) «=* .Igj (control). 

Same as 3 except ai = 2.5xl0 -2 . 

Same as 3a except T g — T g (control) + 10. 

Same as 3b except T g = T g (control) H- 10. 


c 

d 

Figure 8a 
b 
c 

Figure 9 

Figure 10a 
b 
c 
d 
e 
f 

Figure 11 
a 
b 
c 

Figure 12 
a 
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Same as 3a except T g - T g (control) ■ 10. ORIGINAL 
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Same as 3b except T x - 7^ (control) - 10. 

Same as 3b except cloudwater is totally precipitated. 

Same as 4b except cloudwater is totally precipitated. 

Relative humidity field for case where o>/10 is used and cloudwater totally pre- 
cipitated. 

Same as 3 except Q - c 3 *10 2 . 

Same as 3a except 0 C - 9 e (control) + 10 and 7^ — T g (control) 4- 10. 

Same as 3b except 0 e — 0 e (control) + 10. 

Same as 10b except for experiment in which all cloudwater is precipitated. 
Same as 3a except 0 C “ 0 e (control) - 10 and T s *=* T g (control) - 10. 

Same as 3b except 0 e ■» 0 C (control) - 10 and T g — T g (control) - 10. 

Same as lOe except for experiment in which all cloud water is precipitated. 
Control run for model with cloud radiative feedbacks, 
cloudwater with contours at 0, 1, 2, 4, 8 x 10~ 5 kg/kg. 
relative humidity field. 

potential temperature field minus the zonal mean temperature. 

Control run for model with cloud radiative feedbacks. 

infrared heating rates, — , with contours every 4, 2, 1, 0,-1, -2, -4, -8 x 
10“ 5 K 

solar heating rates, with contours every .25, .5, .75, 1 x 10~ 5 K s~ l . 


Figure 13 
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Surface fluxes in watts m~‘ 1 . Positive values increase surface temperature and 
negative values decrease surface temperature. Fyj, 5, SH,E , Frt refer to 
downward thermal radiation, downward solar radiation, sensible heat, 
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evaporation, and upward thermal radiation. Also shown are the surface tem- 
perature T g and 'planetary albedo a. 

Figure 14 Change in temperature after 30 days for the experiment with no cloud radiation 
feedbacks, A0,v, and no solar-cloud radiation feedbacks, M NS . The latter 
experiment includes infrared feedbacks. 

Figure 15 Fixed cloud, FC, and variable cloud, VC, temperature changes due to C0 2 

quadrupling, as a function of time. Both experiments were initialized from the 
end of the control run. 

» 

Figure 16 Differences in cloud field for VC C0 2 experiment minus control. Only plus 
and minus values are indicated because of the large variance of cloud water 
concentrations. In the negative region at the top of the cloud the maximum 


cioud water change is — l.lxlO" 5 kg/ kg. 

Figure 17 Change in planetary albedo for sensitivity experiments minus the control run. 

With the exception of one smalt area for the valuable cloud experiment, the 
albedo change is negative indicating that model atmosphere absorbs more solar 
radiation with increased C0 2 . 


Figure 18 Vertical structure of the potential temperature change (perturbed minus con- 
trol) induced by C0 2 quadrupling for the fixed cloud experiments at 60, 120, 
180, 360 days (FC60, FC120, FC180, FC360) and the variable cloud experi- 
ments at 60, 120, 180 days (VC60, VC 120, VC180).. 
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